Dynamics of a BEC bright soliton in an expulsive potential 
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We theoretically investigate the dynamics of a matter-wave soliton created in a harmonic potential, 
which is attractive in the transverse direction but expulsive in the longitudinal direction. This 
Bose-Einstein-condensate (BEC) bright soliton made of '^Li atoms has been observed in a recent 
experiment (Science 296, 1290 (2002)). We show that the non-polynomial Schrodinger equation, an 
effective one-dimensional equation we derived from the three-dimensional Gross-Pitaevskii equation, 
is able to reproduce the main experimental features of this BEC soliton in an expulsive potential. 



PACS Numbers: 03.75.Fi, 05.45.Yv 
I. INTRODUCTION 

Bright and dark solitons are localized waves that travel 
over large distances without spreading. Bright solitons 
are local maxima in the field while dark solitons are local 
minima. Solitons are ubiquitous: they appear in systems 
as diverse as oceans, shallow waters in narrow channels, 
electric circuits and optical fibers [1]. Recently solitons 
have been experimentally produced in Bose-Einstein con- 
densates (BECs) of alkali- metal atoms [2-5]. Dark soli- 
tons have been obtained with repulsive ®^Rb atoms [2], 
while bright solitons have been observed with attractive 
^Li atoms [3,4] and also in a optical lattice with ®^Rb 
atoms (gap bright solitons) [5]. 

BEC solitons have been theoretically investigated in 
various papers but only in few of them there is a de- 
tailed comparison with experiments [6-8]. For instance, 
by using the Gross-Pitaevskii equation [9], in Ref. [8] we 
have successfully simulated the formation and dynam- 
ics of the train of bright solitons observed in the ex- 
periment of the Rice University [3]. In this paper we 
consider instead the experiment of the Ecole Normale 
Superieure [4], where Khaykovich et al. have produced 
and detected a single BEC bright soliton in a expul- 
sive potential. The stable configurations of this BEC 
bright soliton have been analyzed by Carr and Castin [7] 
by means of an analytical variational method. Here we 
complete their theoretical investigation by studying nu- 
merically the full time-evolution of the bright soliton by 
using the non-polynomial Schrodinger equation (NPSE) 
[10], an effective one-dimensional (ID) equation we have 
recently derived from the 3D Gross-Pitaevskii equation 
[9]. We find that NPSE reproduces remarkably well the 
experimental data, in particular the center of mass dy- 
namics of the BEC cloud and also the time dependence 
of its longitudinal width. 



II. EXPULSIVE POTENTIAL AND NPSE 

In the experiment of Khaykovich et al. [4] the Bose- 
Condensed ''Li atoms are confined in a trap that can be 
modelled by a harmonic potential 
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where m is the atomic mass, uj± — 2t: x 710 Hz is the 
transverse frequency and ujz = 27ri x 78 Hz is the imagi- 
nary longitudinal frequency. This imaginary frequency is 
due to an offset magnetic field that produces the slightly 
expulsive harmonic potential —{m/2)\uJz\'^z'^ in the longi- 
tudinal direction. The main effect of this expulsive term 
is that the center of mass of the BEC accelerates along 
the longitudinal direction [4] . 

The dynamics of a BEC at zero temperature can be de- 
scribed by the time-dependent 3D Gross-Pitaevskii equa- 
tion (3D CPE), given by 
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where ip(r,t) is the macroscopic wave function (order 
parameter) of the Bose-Einstein condensate [9], g — 
ATih^as/m is the interatomic strength, as is the s-wave 
scattering length, and N is the number of condensed 
atoms. 

In the case of strong cylindric radial confinement, 
i.e. when the BEC travsverse energy E±_ of the BEC 
is equal to the transverse harmonic energy Tiui^^ the 
3D GPE reduces to a ID GPE. This result can be 
easily obtained with a variational approach by us- 
ing the following trial wave function ip{x,y,z,t) — 
f{z, t) exp [— (x^ -I- y'^)/{2rf)]/ y/inf, where f{z, t) is the 
longitudinal wave function normalized to one and rj is the 
transverse width. Actually one finds 
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where E± — {fiLo±/2){a\/rj^ -f rj^ /a\) is the transverse 
energy of the Bose condensate. In the ID GPE one has 
rj = a± for which E± = ?iuj±. Instead, choosing a space- 
time dependent ry, one finds [10] that rj satisfies the equa- 
tion 
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V{z,t)=a^ {l + 2asN\f{z,t)\^y/^ 
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By inserting this formula in the previous differential 
equation one gets a nonpolynomial Schrodinger equation 
(NPSE) [10], that reduces to the ID GPE in the weakly- 
interacting Umit a.sNlfl'^ « 1. The NPSE has been 
found to be very accurate in the descripition of BECs un- 
der transverse harmonic confinement and a generic lon- 
gitudinal external potential. The NPSE is accurate with 
both repulsive and attractive scattering length [8,10,11]. 
In our last paper [12] wc have extended the NPSE to in- 
clude also beyond mean- field efi'ects, like the formation 
of a gas of impenetrable bosons in the Tonks-Girardeau 
regime. 



III. BEC BRIGHT SOLITONS 

For simplicity we consider first the case with lOz = 
in Eq. (3). In this case, it is well known that the ID 
GPE (namely the Eq. (3) with r] = a±) with negative 
scattering length (a^ < 0) admits bright soliton solutions 
[1], that set up when the negative inter-atomic energy of 
the BEC compensates the positive kinetic energy such 
that the BEC is self-trapped. Scaling z in units of a± 
and t in units of wj^, with the position 
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from ID GPE one finds the text-book bright soliton 
^{z — vt) = sech ['y{z — vt)] , 



(6) 



where 7 = \as\N/a± and the chemical potential /i is given 
by = 1 — 7^/2, while the velocity v of the bright soliton 
remains arbitrary. Given a value of the scattering length 
Us, the number N of condensed atoms fixes the chemical 
potential ^ of the bright soliton. 

The solitary wave solution of Eq. (6), that we call ID 
bright soliton, exists for any positive value of 7. On the 
other hand, it has been theoretically predicted [10,13] 
and experimentally shown [3,4] that above a critical in- 
teraction strength the BEC bright soliton, that we call 
3D bright soliton, will collapse. To study the properties 
of a 3D bright soliton it is not necessary to solve the full 
3D GPE; in fact, the 3D bright solitons can be very ac- 
curately described by NPSE [8,10]. In particular, from 
the NPSE (namely the Eq. (3) with 77 given by Eq. (4)) 
one finds the 3D bright soliton written in implicit form 
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It is important to observe that this equation is well de- 
fined only for 7$^ < 1/2; at = 1/(27) the transverse 
width ri becomes zero. Moreover, by imposing the nor- 
malization condition, one finds 
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(l-/.)V2 + _^ = 0, (8) 



which admits solutions only for < 7 < 2/3. 

There are thus two kind of collapse: the transverse 
collapse, when the condensate probability density p = 
|$|^ exceeds pc = 1/(27), and the 3D collapse, when 7 
reaches the value 2/3. Note that for a single 3D bright 
soliton at 7c — 2/3 the transverse size of the condensate 
soliton is not yet zero, in fact 7cp < 1/2. Nevertheless, if 
we consider two colliding 3D bright solitons with 7 < 2/3, 
at the impact the total density of the cloud can become 
equal to pc and the system will collapse. 
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FIG. 1. Maximum probability density pm of bright soli- 
tons as a function of the nonlinear strength 7 = \as\N/a±. 
Length z in units Uz = (h/muiz)^^^ , density p in units l/a±, 
where a± = {h/mui±y^^ is the harmonic length of transverse 
confinement. 

In Figure 1 we show the maximum probability density 
Pm of bright solitons as a function of the strength 7 . The 
solid line is the density of a single 3D bright soliton that 
ends at 7 = 2/3. The dot-dashed line is the probability 
density of the ID bright soliton, that coincides with the 
solid line for small values of 7 (in fact, for 7 <C 1 Eq. 
(7) reduces to Eq. (6)). The dashed line is the maxi- 
mum probability density of two colliding equal 3D bright 
solitons, that ends when the curve meets the dotted line, 
which is the density Pc = 1/(27) of the transverse col- 
lapse. Figure 1 shows that two colliding and equal 3D 
bright solitons produce a transverse collapse of the con- 
densate for 7 = 0.472, which corresponds to a single soli- 
ton probability density equal to 0.265. Only for lower 
values of 7 the collapse is avoided during collision and 
interference. 
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IV. BRIGHT SOLITON IN THE EXPULSIVE 
POTENTIAL 

We now consider the case ^ in Eq. (3) and try to 
simulate the experiment of the Ecole Normale Superieure 
[4] with ^Li alkah-metal vapors. In that experiment a 
BEC with positive scattering length (a^ > 0) has been 
condensed into a cigar-shaped harmonic trap. The scat- 
tering length has been then tuned to a small and negative 
value {as < 0) via a Feshbach resonance. The resulting 
condensate, projected onto the expulsive harmonic po- 
tential of Eq. (1), has been observed to propagate over a 
distance of 1 mm. 




FIG. 2. Density of the '^Li BEC in the expulsive poten- 
tial obtained by solving the NPSE. The BEC cloud propa- 
gates over 1 mm. Case with as = (ideal gas). There are 
A'' = 4 X 10^ atoms. Six frames from bottom to top: t = 2 ms, 
t — 3 ms, t — 4 ms, t — 5 ms, t — 6 ms, t = 7 ms. External 
harmonic potential given by Eq. (1). Red color corresponds 
to highest density. 

Following the paper of Khaykovich et al. [4] , we choose 
as initial condition for a fixed the ground-state of the 
stationary Eq. (3) with a fully confining harmonic poten- 
tial: wj^ = 27r X 710 Hz and = 27r x 50 Hz. Then the 
time-dependent wave function is obtained by the numeri- 
cal integration of the Eq. (3) with the expulsive potential: 
= 27r x 710 Hz but = 27ri x 78 Hz. Note that, as in 
the experiment, the initial position of the BEC cloud is 
shifted of 50 /im on the left with respect to the maximum 
of the expulsive potential. A critical point is the choice of 
the number N of atoms. To avoid the collapse we choose 



A'' = 4 X 10'^, a value compatible with the experimental 
data and suggested by the variational theory [13]. The 
NPSE is solved by using the finite-difference numerical 
algorithm described in [14]. 




FIG. 3. Density of the Li BEC in the expulsive potential 
obtained by solving the NPSE. The BEC cloud propagates 
over 1 mm. Case with as = —0.21 nm ("bright soliton"). 
There are A'^ = 4 x 10'^ atoms. Six frames from bottom to 
top: t = 2 ms, t = 3 ms, t = 4 ms, t = 5 ms, t — 6 ms, t — 7 
ms. External harmonic potential given by Eq. (1). Red color 
corresponds to highest density. 

In Fig. 2 we plot six frames of the density of the BEC 
with Us — 0. This is the case of an ideal gas and the 
figure shows that the center of of the bosonic 

cloud follows the law Zcm{t) = exp (Iw^li). In addi- 
tion, the longitudinal width of the cloud grows due to 
the dispersive kinetic term. In Fig. 3 we plot instead 
the time-dependent density of the BEC with as = —0.21 
nm. Also here the center of mass Zcm follows the expo- 
nential grow in time but the longitudinal width does not 
show an appreciable enlargement: the matter wave trav- 
els without spreading. It is important to stress that Fig. 
2 and Fig. 3 are in remarkable good agreement with the 
absorption images shown in Fig. 3 of the experimental 
paper [4] with the same values of the parameters. 

In their experiment, Khaykovich et al. [4] have mea- 
sured the root mean square size a of the longitudinal 
width versus the propagation time for three values of a^: 
as — 0, Qs = —0.11 nm, and a^. = —0.21 nm. Figure 
4 shows their experimental data and our numerical re- 
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suits obtained with the NPSE. For ag = the attractive 
interaction between atoms is zero and the expansion of 
the cloud is governed by the kinetic energy and the neg- 
ative curvature of the longitudinal axial potential. For 
Us = —0.11 nm the size a of the axial width is consistently 
below that of a non-interacting gas but the attractive in- 
teraction is not strong enough to avoid an appreciable 
enlargement of the cloud. 
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FIG. 4. Root mean square size a of the longitudinal width 
of the BEG as a function of the propagation time t. The 
filled circles are the experimental data taken from Ref. [4]. 
The dashed line is the ideal gas (a^ = 0) curve. The solid line 
is obtained from the numerical solution of the NPSE. 

For fls = —0.21 nm the enlargement is further reduced 
but our numerical results suggest that a is not truly 
constant. This is not surprising because a truly shape- 
invariant solitary wave is expected only in the absence of 
a longitudinal potential and with an appropriate initial 
condition. Nevertheless Fig. 4 shows that our results are 
fully compatible with the experimental data, which have 
the resolution limit Act = 9 /xm of the imaging system. 

As previously discussed, in the calculations we have 
choosen = 4 x 10^ to avoid the collapse of the BEC 
with = —0.21 nm. In the experiments the collapse 
actually implies that a fraction of atoms is expelled from 
the Bose condensate via three-body recombination. This 
effect can be phenomenologically modelled by including 
a dissipative term in the GPE or in the NPSE, but a 
more rigorous treatment requires the solution of coupled 
time-dependent equations for the Bose condensate and 
the thermal cloud. 

Finally, we note that a sudden shift of the inter-atomic 
strength \ai\N ja^^ to a very large value (for instance via 
a Feshbach resonance) produces a set of bright solitons 



(soliton train), which can travel in the expulsive poten- 
tial. A detailed analysis of the formation of these multi- 
soliton configurations can be found in our recent papers 
[8,15]. 



V. CONCLUSIONS 

We have shown that the ID non-polynomial 
Schrodinger equation (NPSE) is able to accurately de- 
scribe the time evolution of an attractive Bose-condensed 
cloud of ^Li atoms in an expulsive harmonic poten- 
tial. Contrary to the ID Gross-Pitaevskii equations, the 
NPSE takes also into account the space-time variations 
of the transverse width of the Bosc-Einstcin condensate 
and its numerical integration is much faster than the 
numerical solution of the 3D Gross-Pitaevskii equation. 
The theoretical results of the NPSE have been compared 
with the experimental data of a recent experiment per- 
formed at the Ecole Normale Superieure of Paris. Both 
experiment and theory show that, when the inter-atomic 
strength is sufficiently strong, a localized matter wave is 
obtained, which travels for a large distance without an 
appreciable spreading. Our numerical calculations sug- 
gest that this solitary wave is not fully shape-invariant 
but the experimental resolution does not allow one to 
clarify this point. It is important to stress that the limita- 
tions in the current resolution of the experimental imag- 
ing system can be overcome by creating a bright soliton 
with a larger longitudinal width. This can be achieved 
by using a smaller attractive inter-atomic strength in ad- 
dition with a weaker expulsive potential. 

In conclusion, we observe that to make clear the 
particle-like nature of a soliton it is essential to inves- 
tigate its interaction with another soliton. The interfer- 
ence of BEC bright solitons is an important subject which 
has to be analyzed in detail, having also implications for 
atom interferometry [16]. 
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